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04 ' We study the prospects for observing H2 and HD emission during the assembly 

^ I of primordial molecular cloud cores. The primordial molecular cloud cores, which 

resemble those at the present epoch, can emerge around 1 + z ~ 20 according to 
recent numerical simulations. A core typically contracts to form the first generation 
of stars and the contracting core emits H2 and HD line radiation. These lines show a 
double peak feature. The higher peak is the H2 line of the J = 2 — (v=0) rotational 
transition, and the lower peak is the HD line of the J = 4 — 3 (v=0) rotational 
transition. The ratio of the peaks is about 20, this value characterising the emission 
from primordial galaxies. The expected emission flux at the redshift of 1 + z ~ 20 
(e.g. rim — 0.3 and flA = 0.7), in the J = 2 — (v=0) hue of H2 occurs at a rate 
~ 2 X lO"'^ Jy, and in the J = 4 - 3 (v=0) hne of HD at a rate ~ 8 x 10"^ Jy 
P ■ The former has a frequency of 5.33179x10" Hz and the latter is at 5.33388 x 10"Hz, 

respectively. Since the frequency resolution of ALMA is about 40 kHz, the double 
^ • peak is resolvable. While an individual object is not observable even by ALMA, the 

^ ' expected assembly of primordial star clusters on subgalactic scales can result in fluxes 

at the 2000-50 /zJy level. These are marginally observable. The first peak of H2 is 
produced when the core gas cools due to HD cooling, while the second peak of HD 
rS I occurs because the medium maintains thermal balance by H2 cooling which must be 

^ . enhanced by three-body reactions to form H2 itself. 

Key words: cosmology: observations — galaxies: formation — ISM: molecules — 
submillimeter 



1 INTRODUCTION 

Observation of the first generation of stars presents one of 
the most exciting challenges in astrophysics and cosmology. 
Hydrogen molecules (H2S) play an important role as a cool- 
ing agent of the gravitationally contracting primordial gas 
(Saslaw & Zipoy 1967) in the process of primordial star for- 
mation. It has been argued that H2 is an effective coolant for 
the formation of globular clusters (Peebles & Dicke 1969), 
if such objects precede galaxy formation. The contraction 
of primordial gas was also studied in the pioneering work 
of Matsuda, Sato, & Takeda (1969). Thus, we expect the 
detection of H2 hne emission from primordial star-forming 
regions (Shchekinov 1991; Kamaya & Silk 2002; Ripamonti 
et al. 2002). Observational feedback from such Population 
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III objects was examined in Carr, Bond, & Arnett (1984). 
In the context of the CDM (Cold Dark Matter) scenario 
for cosmic structure formation, Couchman & Rees (1986) 
suggested that the feedback from the first structures is not 
negligible for e.g. the reionisation of the Universe and the 
Jeans mass at this epoch (e.g. Haiman, Thoul, & Loeb 1996; 
Gnedin & Ostriker 1997; Ferrara 1998; Susa & Umemura 
2000). Recent progress on the role of primordial stars for- 
mation in structure formation is reviewed in Nishi et al. 
(1998) and its future strategy is discussed in Silk (2000). 

More recently, the first structures in the Universe have 
been studied by means of very high resolution numerical 
simulations (Abel, Bryan, & Norman 2000). The numerical 
resolution is sufficient to study the formation of the first 
generation of molecular clouds. According to their results, a 
molecular cloud emerges with a mass of ~ 10^ solar masses 
as a result of the merging of small clumps which trace the 
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initial perturbations for cosmic structure formation. Due to 
tiie cooling of H2, a small and cold prestellar object appears 
inside the primordial molecular cloud. It resembles the cores 
of molecular clouds at the present epoch. These numerical 
results are consistent with other numerical simulations by 
Bromm, Coppi, & Larson (1999; 2001). All the simulations 
predict that the primordial molecular clouds and their cores 
appear at an epoch of 1 + z ~ 20 (2: is redshift). 

According to these results, the first generation of young 
stellar objects has a mass of ^ 200 solar masses, density 
of ~ 10^ cm~^, and temperature of ~ 200 K. We stress 
that these are cloud cores, and not stars. Inside the cores, 
a very dense structure appears. We call this a kernel for 
clarity of presentation (Kamaya & Silk 2001). Its density 
increases to a value as high as ~ 10* cm~"^ where three-body 
reactions for H2 formation occurs (Palla, Salpeter, & Stabler 
1983, Omukai 2000). Further evolution is by fragmentation 
(Palla, Salpeter, & Stabler 1983) and/or collapse (Omukai 
& Nishi 1998). In the previous paper, we considered the 
collapsing kernel and the possibility for its observation by 
future facilities. H2 line emission tracing the temperature 
structure of the kernel is potentially detected by ASTRO- 
and ALMA|^ if the kernels are collectively assembled, 
as might be expected in a starburst. We also comment here 
that if there is dust grain in starforming region, the situation 
changes (e.g. Hirashita, Hunt, Ferrara 2002). For example, 
the three body reaction dominates over the grain surface 
formation if dust grain is present type and only when gas 
density is above 10^^ cm~^. 

Now in the low temperature regime of the typical core, 
HD is also an important coolant (e.g. Shchekinov 1986). In- 
deed, Uehara & Inutsuka (1998) show numerically that the 
core collapses at a relatively low temperature if HD exists. 
According to them, the primordial filament can reach a tem- 
perature of 50 K. Thus, the effect of HD should be examined. 
Indeed as long as HD is a dominant coolant, it is useful to 
consider the observational possibility of HD emission from 
primordial molecular cloud cores, as we do below. 

Furthermore, because of the large electric dipole mo- 
ment, LiH is also potentially an important coolant (Lepp & 
ShuU 1984). Despite the very low abundance of LiH, which 
means that it is never the dominant coolant, as shown in 
Lepp & ShuU, if its lines are optically thin, it can still be 
important. We also comment on line cooling by LiH in this 
paper. 

In addition to a contracting primordial cloud, another 
molecular emission line is discussed by Ciardi & Ferrara 
(2001) who predict the detection of mid-infrared emission 
from primordial supernova shell by NGST|^ According to 
Ciardi & Ferrara, the shell can cool and radiate in rotational- 
vibrational emission lines, and then mid-infrared emission is 
expected. In particular, the shell can also emit the possible 
ground rotational line of H2 (J = 2 — and t; = 0). In our 
first paper, we also predict the primordial emission of the 
ground rotational line of H2. If and when primordial emis- 
sion of the ground rotational transition is detected, it will 
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be important to judge which is its origin. In this paper, we 
present a very simple answer to this question. 

In §2, we formulate the cooling function for HD and 
LiH. For H2, this was presented in a previous paper (Kamaya 
& Silk 2001). We summarise the cooling function of H2 in 
the Appendix. In §3, we present a model structure for the 
primordial molecular cloud cores. In §4, we estimate that the 
luminosity from the cores is primarily due to the emission of 
H2, with a secondary role played by HD, and that LiH is not 
important. In §5, the observational feasibility of detection 
is reviewed (differently from the discussion in our previous 
paper). The theory underlying the observational predictions 
is presented in §6, and our paper is summarised in the final 
section. 



2 DESCRIPTION OF MOLECULAR EMISSION 

Primordial molecular cloud cores are expected to be 
low temperature objects (Abel et al. 2000), where HD 
(Shchekinov 1986) and LiH (Lepp & ShuU 1984) line cool- 
ing are frequently considered to be important. In this sec- 
tion, first of all, we present a formalism for the molecular 
line emission to estimate the luminosity from the cores. The 
lower rotational transition emission of H2 is also important. 
Since a detailed description has been presented in the pre- 
vious paper (Kamaya & Silk 2001), only a brief summary is 
given in the Appendix. 



2.1 HD 

HD has a weak electric dipole moment, since the proton in 
the molecule is more mobile than the deuteron; the elec- 
tron then does not follow exactly the motion of the positive 
charge, producing a dipole (Combes & Pfenniger 1997). This 
moment has been measured in the ground vibrational state 
from the intensity of the pure rotational spectrum by Tre- 
fler & Gush (1968). The measured value is about 5.85 10~^ 
Debye (1 Debye = 10~^* in cgs unit). Since the first rota- 
tional level is about 128 K, the corresponding wavelength 
is 112 /.im. Since the temperature of the primordial molec- 
ular cloud cores (Abel et a. 2000) is about 200 K, HD has 
sufficient potential to be the dominant coolant (Shchekinov 
1986). In this paper, we set the D abundance to be 4 x 10~^, 
keeping consistency for our structural model of cores with 
their thermal conditions predicted by Uehara & Inutsuka 
(2000), which is the most recent and detailed analysis for a 
primordial cloud undergoing HD cooling. 

For the cooling function, we adopt the formalism of Hol- 
lenbach & McKee (1979). According to this. 



AnD.thin = nXWHDX 



4(fcT)2ylo 



nEo{l + (ncr/n) + 1.5(ncr/n)0-5) 



erg cm 
(1) 
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where k is the Boltzmann constant, n is the gas number den- 
sity. Tier is the critical density, nno is HD number density, 
Ao is 5.12 X 10"® s"^ (AbgraU et al. 1992), and Eq = 64fc 
erg for HD. Determining ricr exactly is a complex calcula- 
tion, while for our purpose it is sufficient to consider ncr 
for the dominant cooling transition if we are interested in 
the dominant process. Hence, we estimate ricr via the for- 
malism of HoUenbach & McKee (1979) as ricr = 7.7 x 10^ 



Observing the double emission Line feature from primordial cloud cores 3 



cm~"^ (Tg/lOOO.O)""'' where Tg is the thermal temperature 
of gas. For reader's convenience, we would like to note that 
our AnD.thin has changed from L of (6.23) of Hollenbach & 
McKee (1979) by the amount of n x nnD, while both are 
always consistent. 

However, the above formula is inadequate if the temper- 
ature is low since the population of the rotationally excited 
level is small. We consider the cooling function of HD in 
the optically thin and low temperature regime below 255 
K (Galli & Palla 1998) if the density is below the critical 
density for a given temperature; 

Ahd;„^o = 27ioSioexp(-£;io/feTg)-|-(5/3)72iS2iexp(-£;2i//cTg) 

(2) 

Here, Eu) = 128k erg, E21 = 255k erg , 710 = 4.410"^^ + 
S.eiO^^^Tg" '^^ cm^ s~\ and 721 = 4.110"^^ + 2.110-"Tg°-^2 
cm^ s~^. If the gas density is above the critical density, we 
simply estimate the cooling function by two methods. The 
first is for the temperature range from 128 K to 255 K. In this 
range, we simply use Eq.(l) multiplied by exp(— 256K/Tg). 
This is because the population of J = 2 to that of J = 1 
is reduced by the factor of exp(256K/rg). The second is for 
if the temperature is below 128 K. In such low tempera^ 
ture case, as simple modification as for the first method for 
correction is not applicable. Fortunately, the two-level sys- 
tem is a good approximation. Then, we estimate the cooling 
function from the spontaneous de-excitation rate (A) of the 
most probable rotational level of J = 1 and the coUisional 
de-excitation rate. The de-excitation rate is estimated from 
the coUisional cross section of dc-cxcitation divided by a 
sound speed. The coUisional cross section of de-excitation 
is typically 2.0 x 10"^^ cm^ for HD (Schaefer 1990). Our 
approximation breaks at the outer edge of the core, where 
we smoothly extrapolate keeping the thermal balance. For- 
tunately, this correction docs not alter the main conclusion 
and the dominant cooling level of J. 

On the other hand, to find the cooling function in the 
optically thick regime, we always need the Einstein A coef- 
ficient. We calculate it approximately in the standard way. 
The coefficient is 



J 



2J + 1 



(3) 



where Vjij is the frequency of the J' J transition, h 
the Planck constant, c the light speed, D the electric dipole 
moment, and J the rotational quantum number. The results 
via Eq.(3) describe the exact values of Abgrall et al. (1982) 
for HD very well. 

Once Jd(r(r)) which is the dominant cooling level of 
J is determined at each position of the core, we can mod- 
ify the optically thin cooling function to the optically thick 
cooling function by multiplication of the escape probability, 
e. To obtain the cooling function for the thick case, hence 
we define: 



1 -exp(r.,,,) 



and 



9j' 

gj 



n.j' \ Rj 
nj J 6v 



(4) 



(5) 



and g,j the statistical weight of 2 J -|- 1. The velocity disper- 
sion corresponds to Doppler broadening, and is estimated 
to be given by the sound speed. The inferred optically thick 
cooling function is then 



AhD, thick = AHD.thin X ivjij- 



(6) 



The procedure of Jeans length shielding in (5) is useful for 
a simple analytical analysis (Low & Lynden-Bell 1976; Silk 
1977). 

2.2 LiH 

The LiH molecule has a much larger dipole moment of 5.9 
Debye and the first rotational level is only at 21 K. Thus, 
although its abundance is very small (~ 10"^°; we adopt 
this value), there is a plenty of possibility for LiH to be 
an important coolant (Lepp & ShuU 1984). For LiH, fortu- 
nately, the formalism of Hollenbach & McKee (1979) is a 
very good approximation. Then, we use Eq.(l) with param- 
eters for LiH in the thin regime, and the escape probability 
correction is applied if the optical depth of the dominant 
line emission is above unity. The adopted parameters are 
Ao = 0.0113 s"^ and Eo/k = 21.0 K, and "LiH instead of 
riHD. For the corresponding de-excitation rates, we use the 
results of detailed balance analysis (Bougleux & Galli 1997). 
The expected coUisional do-oxcitation cross section is found 
to be approximately 5.6 x 10"^'' cm"^, then the de-excitation 
rate becomes 2 x 10"^° cm^ s"^ at Tg = 3000 K. Since 
the de-excitation cross sections tend to be nearly constant 
at low energies for collisions with neutral particles, the ap- 
proximate de-excitation rate at lower temperature can be 
obtained multipUed by the factor of (rg/3000 K)" ^ If the 
comparison with the case of collision of He are interested, 
the reduced factor of 3*^^ is adopted to account for the dif- 
ference in the reduced mass (Bougleux & Galli 1997). Their 
appendix B is useful for further details. If only the optically 
thin cooling function is needed, appendix A. 3 of Galli & 
Palla (1998) is useful as long as the gas density is below any 
critical density. 



3 STRUCTURAL MODEL OF PRIMORDIAL 
MOLECULAR CLOUD CORE 

According to recent numerical simulations, molecular cloud 
cores appear prior to the formation of population III objects. 
Gores contain a dense and cool kernel. In our previous paper, 
we consider how much the kernel emits H2 line luminosity. In 
the current paper, we consider the case of the cores. Since the 
cores have lower temperature than the kernels (e.g. Uehara 
& Irmtsuka 2000), HD and LiH nmst be cousidered. 

The emission properties are determined by the density 
and temperature structure of the cores. Fortunately, a rea- 
sonably simple model is possible according to Uehara & In- 
utsuka (2000) and Omukai (2000). We find a fitting formula 
for the distribution of H2 and temperature which is found by 
trial and error to approximately reach thermal equilibrium 
(see the next section). For f2{r) (solid line in Fig.l); 



/2(r) = 0.0001-F0.495X 



exp 



n(r) ^ 
lOll Ocm-3 J 



exp I 



n(r) > 
lOll-Ocm-3 ; 



where Rj is the Jeans length, 5v is the velocity dispersion. 



e^p( ioi"ol-0 +e^p(ToT^ 



n(r) 



(8) 
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where n(r) = 10.0*^cm^'^(r/0.01 pc)^'^ '^ (the maximum of 
/2(r) is set to be 0.5 by definition). This describes approxi- 
mately the effect of the three body reaction to form H2. For 
T(r) (dashed Une in Fig.l); 



T(r) = 50 K 



n(r) 



1040cm-3 



(9) 



We hypothesise that all D and Li are in molecular form 
above a density of 10'* cm~^. This is partially supported by 
Uehara & Inutsuka (2000). 

To examine the total emission energy, we need to de- 
termine the mass distribution around the centre of the core, 
where the first star emerges. Also, we assume a spherical con- 
figuration for the mass distribution. According to Omukai & 
Nishi (1998), a high accretion rate is realized if a similarity 
collapse occurs with the adiabatic heat ratio of 1.1 (e.g. Suto 
& Silk 1988). The density distribution is described as 

Slnp(r) _ -2 



dlnr 



(8) 



2-7 

Here, r is the radial distance from the centre, p{r) is the 
mass density of atomic H, H2 and He, and 7 is the spe- 
cific heat ratio. We set the mass-density distribution of a 
protostellax-core with 7 = 1.1 as p(r) = po{r /ro)~^'^ where 
po is 2.0x10"^° g cm~^ and ro is 0.63 pc. These values are 
appropriate for fitting a typical protostellar core of Abel et 
al. (2000). 



4 MOLECULE EMISSION LUMINOSITY 

We are interested in the dominant rotational emission, then 
we calculate 

/ 47rr2A,,thick(r)Jd(r(r))dr ^ ^ 

Jmax = \. . . ■ (10) 



Here, for HD (z=HD) 

Jd(T(r)) = Jd,HD(T(r))= fjx ^^''^ 



2 64.0 K 



and for LiH (i=LiH) 



Jd(r(r)) = Jd,LiH(T(r)) = ( I X ^('■^ 



and for H2 (i=H2) 

Jd(r(r)) = Jd,H,(T(r)) = 



2 21.0 K 



T{r) 



85.0 K 



(11) 



(12) 



(13) 



Here, all Jd{T{r)) mean the dominantly contributing J-level 
to the statistical weight at temperature of T{r) (Silk 1983). 
Using our fitting formula of /2(r) and T{r), we find each 
Jmax is about 4.0 for HD, 17.0 for LiH, and 2.3 for H2, 
respectively. Then, to estimate the dominant line luminosity 
of Lthick of Eq.(6), we use J = 4 for HD, J = 17 for LiH, 
and J = 2 for H2, respectively. The exceptional treatment 
(but reviewed in Kamaya & Silk 2001) for H2 is given in 
the Appendix, according to which we can discuss J = 2 — 
transition of H2 independent of any other transition. 

The estimated total luminosity for a single core 
(/core ^Trr^Athickdr), Lsingie, is 4.2 x 10^® erg s"^ Each com- 
ponent is 4.1 X 10^^ erg s"^ for H2, 2.1 x 10^"* erg s"^ for HD, 



and 6.9 x 10'^" erg s^^ for LiH. Here, first of all, we find that 
the contribution due to LiH is very small. This contradicts 
the conclusion of Lepp & ShuU (1984). Fortunately, the rea- 
son is simple. Although Lepp & Shull regarded all LiH lines 
as optically thin, they are optically thick for our case in the 
region where LiH would be important £is suggested by Lepp 
& Shull. 

The line broadening is estimated to be ~ 5vr> = 
(2fcT//i(r)mH)°'^ in the dimension of velocity. Here, p,(r) is 
mean molecular weight at each position. Adopting this, the 
luminosity per Hz (T=1000 K is assumed) is as follows; H2 
rotation emission of J = 2-0 (1.1x10^^ Hz; 2.8x10"^ mm) 
is 4.0 xlO^'' erg Hz"^, HD rotation emission of J = 4 - 3 
(1.1x10^-^ Hz; 2.8x10"^ mm) is 2.0 xlO^® ergHz"\ and LiH 
rotation emission of J = 17-16 (7.4x10^^ Hz; 4.0x10"^ 
mm) is 9.8 xlO^^ erg Hz~^. Here, we consider the dominant 
Jmax to Jmax — 1 transitions for HD and LiH, and Jmax to 
Jmax — 2 transition for H2. 

The total cooling rate at each position is summarised in 
figure 2. Below 10* cm~'^, HD is the dominant coolant, while 
H2 is dominant above this density. This confirms the results 
of Lepp and Shull (1984). We also check if our model of the 
molecular cloud core is reasonable or not. To do it, we esti- 
mate the heating rate. In a contracting cloud without dust, 
the compressional heating is generally dominant. When the 
three-body reactions for H2 formation occur, chemical heat- 
ing is also important. We consider these two heating mech- 
anisms. For compressional heating, we estimate c^/tg (e.g. 
Omukai 2000). The specific free energy for the three-body 
reaction is 4.48 eV. Our result is displayed in figure 3. The 
ratio of cooling and heating is given in figure 4. As clearly 
shown in this figure, the deviation from the thermal balance 
is within a factor of three. Thus, our model structure for the 
core is consistent with thermal balance between molecular 
line cooling and the expected heating over all of the density 
range considered in this paper. 



5 OBSERVATIONAL FEASIBILITY OF 
DOUBLE PEAK EMISSION 

According to our first paper (Kamaya & Silk 2001), some 
bright emission lines from assembly of primordial young stel- 
lar objects usually have sub-mjy flux as long as the redshift 
of z is about 10-40. Then, when we are interested in the 
same range of the redshift, it is sufficient to discuss only a 
typical case of 2 = 19. The expected emission of the same 
emission lines from different redshifts can also have the simi- 
lar flux level. This realises because the flux per frequency has 
apparently positive effect of the redshift (e.g. Eq.(6) of Cia- 
rdi & Ferrara 2001) and the sound speed at the primordial 
starforming region should also be redshifted (Toleman 1966, 
Kamaya & Silk 2001). Furthermore, our estimate is reason- 
able if we do not consider smaller z than 6 at which the 
reionisation of the Universe occurs and the adopted assump- 
tions for analysis break. 

In the current paper, we need the redshifted observa^ 
tional frequency and wavelength. For the three lines, we 
obtain 0.53x10^^ Hz; 5.6x10"^ mm for J = 2 - of H2, 
0.53x10^^ Hz; 5.6x10'^ mm for J = 4 - 3 of HD, and 
3.7x10^^ Hz; 8.1x10"^ mm for J = 17 - 16 of LiH. The 
adopted 1 -|- ^; is 20 according to the results of recent nu- 
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merical simulations. Obviously, the predicted frequency is 
located in the range of ALMA. Hence, we discuss the ob- 
servational possibility of detection by ALMA in the current 
paper. ALMA is a ground-based radio interferometric fa- 
cility, and will consist of many 12-m antennas. A detailed 
recent review is found in Takeuchi et al. (2001). According 
to this, the 5a sensitivities at 350 fxm, 450, 650, 850, 1.3 
mm, 3.0 mm are expected to be 390, 220, 120, 16, 7.5, 4.6 
pjy, respectively (8-GHz bandwidth). 

The most prominent feature is predicted to be a double 
peak of HD and H2 emission. A schematic view is presented 
in figure 5, but it is depicted at the coordinate of the core 
(i.e. they are not rcdshiftcd). The most obvious feature is 
the difference of the peaks of each line intensity. The ratio 
of H2 to HD is about 20. Detection of this double peak fear 
ture would confirms the presence of primordial molecules in 
forming galaxies. The difference in frequency of both the 
molecules is about 10* Hz. Each of the precise values is 
5.3317 10" Hz for H2 (J=2-0) and 5.3338 lO" Hz for HD 
(J=4-3), respectively. The frequency difference is resolved 
sufiiciently by ALMA since it has frequency resolution of 
about 4 X 10* Hz. 

To estimate the observed flux, we determine the dis- 
tance to the object. Then, we calculate it numerically from 
the following standard formula: 

^''":ho/ (nA + nM(i + 2)-^-»)"--' 

where Ho is the Hubble parameter taken to be 75 km sec~^ 
Mpc~^, Oa is the cosmological constant parameter, and Qm 
is the density parameter. We consider the case JlA + fiM ~ 1, 
since our discussion bases on the numerical results (e.g. Abel 
et al. 2000). Adopting D19, we obtain the observed fluxes 
of each of the lines. The results are summarised in table 
1. Although the line- broadening is estimated from voSvo/c, 
we also correct it for the redshift effect {vo is the central 
frequency). For each of the parameter sets of (f^Aj^M), we 
obtain D19 = 0.62x10^° pc (V^a = 0, = 1), Dig = l.OOx 
10^° pc (Oa = 0.7, = 0.3), and Dig = 4.58 x 10^° pc 
{Qa = 0.9, Qm ~ 0.1). According to Table 1, the rotational 
line fluxes of J = 2-0 (v=0) for H2 are 0.16 fiJy; 0.0083 ^iJy 
for J = 4-3 (v=0) for HD; and 0.000001 /iJy for J = 17-16 
(v=0) for LiH if Qm + = 1. Thus, we conclude that a 
single core is not easily observable even by ALMA. 

However we note that if the cores are collectively assem- 
bled on a sub-galactic scale (Shchekinov 1991), the agglom- 
eration can be detected by ALMA (Kamaya & Silk 2001). 
We shall estimate the number of cores in a primordial galaxy. 
Firstly, we must consider the lifetime of a core able to show 
double peak emission. The required density and temperature 
distribution for such a core is possible when the accretion 
rate of the contracting gas is about 0.01 Mq year~^ (Omukai 
& Nishi 1998; Kamaya & Silk 2001). Then, the Ufe-time of 
the core with double peak emission is about 10* years as 
long as a massive star of ~ IOOMq forms inside the core. 
A W^Mq cloud (Abel et al. 2000) is expected to form 1000 
such cores at a plausible efficiency of 0.1 in mass during its 
life-time. Taking its lifetime to be ~ 10® years as the free- 
fall time of a W^Mq cloud, we find that its luminosity can 
reach 100 Lgingie- 

Next, we consider an entire primordial galaxy with 
10^^ Mq. It may form 10^ supernovae over its entire lifetime. 



We assume that it makes 10^ primordial massive stars, since 
this number of massive stars gives enrichment to roughly 1 
percent of the solar metallicity. If the burst of formation of 
primordial molecular cloud cores occurs in the central region 
of the primordial galaxy (we postulate 1 kpc as the size of the 
core- forming region), then there are 10* such giant molecu- 
lar clouds with 10® Mq in the core-forming region, as long 
as the massive star forms with a high accretion rate of 0.01 
Mq year~^ (Kamaya & Silk 2001). During this phase, the 
cumulative luminosity would be 10*' Lsingic- However, the dy- 
namical time-scale in the core-forming region may be ^ 10^ 
years (e.g. the duration of the starburst). Then, we obtain 
10* Lsingic for the luminosity of the primordial galaxy under- 
going its first star formation burst, conservatively assuming 
^ 10® years as the life-time of giant molecular clouds. Fi- 
nally, for molecular-line emitting proto-galaxies, we obtain 
1.6 mJy for J = 2 - of H2, 0.83 mjy for J = 4 - 3 of 
HD and 0.00001 mJy for J = 17 - 16 of LiH (f^M = 0.3 
and Qa ~ 0.7). The estimated flux levels of H2 and HD are 
consistent with Shchekinov (1991), in which only the low- 
est transition lines for both the molecules were discussed. It 
may also be better to say that these values are optimistic 
values. In more realistic conditions, furthermore, our simple 
time-dependent summation scheme might break down. 

According to the current status of the instrumentation 
for ALMA, 80-890 GHz is the allowed detectable range. It is 
feasible for our prediction of the redshifted emission. Unfor- 
tunately, however, the transmission is bad for the predicted 
feature at 530 GHz because of the atmosphere of the Earth. 
Then, it may be necessary to detect H2 and HD emission 
from a protogalaxy forming larger than 1 + z — 20. Since 
Hm < 1, which seems to be reasonable if S7a = 0.7, the 
formation of the cores is prompted. Then, we can expect 
to detect emission at a larger frequency than 530 GHz, and 
detection should be feasible of the H2 and HD double emis- 
sion from the assembly of the primordial molecular cores. Of 
course, we also expect the emission later than 1 + z — 20, 
since the Universe has been reionised at 2: ~ 6. The lat- 
ter case is also favourable for ALMA because of the better 
sensitivity. 

In the previous paper, we advocated a deep blank field 
survey. Then, in the current paper, we propose another ob- 
servational strategy to detect the double emission feature 
from the primordial molecular cloud cores. Measurement of 
the number counts of submillimeter sources is planned as 
one key galaxy evolution project for ALMA (Takeuchi et 
al. 2001). According to the predictions, many submillimeter 
sources are expected below sub-mjy levels. Once the survey 
is operating, we can utilise the submillimeter source count 
data. Firstly, we pick up the faint sources around QSOs 
since QSOs are expected at high density peaks in the usual 
hierarchical model of cosmic structure formation from pri- 
mordial Gaussian-distributed density fiuctuations. Primor- 
dial galaxies are concentrated around QSOs. Secondly, we 
roughly cheek the spectra. Primordial galaxies do not have 
dust, while evolved systems have dust. This means the con- 
tinuum level of the flux at submillimeter wavelengths is sig- 
nificantly different between primordial and evolved galaxies. 
Unfortunately, it seems that the faint sources are observed 
in the low frequency resolution mode according to Takeuchi 
et al. (2001). But, one can expect to find unresolved double 
peaks as a bump in the observed spectral energy distribu- 
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tion, which is a different feature from the blackbody emission 
by dust. To find the bump due to the double emission lines 
should not be a difficult task. After this simple selection, 
we re-observe the candidates for primordial submillimeter 
sources in the high resolution mode in frequency. Finally, 
it is expected that the primordial cores can be discovered 
showing the double peak spectral features due to H2 and 
HD, in which the ratio of the two peaks, which is about 50, 
confirms the hypothesis of primordial emission. 

Finally, we stress that we arc able to recognise the differ- 
ence between H2 emission from the assembly of primordial 
molecular clouds and that supernova shells in primordial 
starforming regions (Ciardi & Ferrara 2001). According to 
figure 1 of Ciardi & Ferrara, the emission of the ground ro- 
tational line of H2 {v = 0), which is never the dominant 
cooling line because of the high temperature of the shell, 
can have a flux level of sub-mjy to mjy at z ^ 10. This flux 
level is obtained via adoption of the sound velocity of 10 km 
sec~^ in the shell gas. Thus, both our predictions and those 
of Ciardi & Ferrara yield similar flux levels for the same 
molecular line of H2. The difference between the two pre- 
dictions is the HD line to be associated with the H2 line in 
our case. To conclude, if and when the double peak emission 
is found, it would strongly indicate the primordial emission 
from the assembly of the first starforming molecular clouds. 



6 THEORY OF DOUBLE PEAK EMISSION 

How is the double-peak emission produced? We describe a 
theory for double-peak emission from primordial molecular 
cloud cores. In the entire region of the core, the temperature 
is significantly lower than 515 K which corresponds to the 
transition energy of H2 (J=2-0). Such low temperatures are 
possible only when HD cooling is efficient. Then, the effect 
of the large volume of the low temperature region lets the 
first and larger peak be a rotational emission line of H2 
(J=2-0). This is possible because of the thermalisation of the 
gas as long as the gas temperature is above 100 K (i.e. ~ 
2^ X 85 X 2/7; J=2), since the rotational level of J = 2 is 
excited by chance. 

The situation of the second peak is a little more compli- 
cated. We define a characteristic density for the three-body 
reaction to set in as being nthree ^ 10® cm~^. Around the 
region of this density, the cooling rates of HD and H2 are 
comparable. Once it is confirmed that the temperature is de- 
termined by the balance of the cooling of HD and H2 with 
adiabatic heating, we find that a temperature of ~ 300 K 
at a density of Wthroc is realized. In other words, if the cool- 
ing of H2 were to be neglected, the temperature would be 
lower than 300 K since the corresponding heating was not 
allowed there. By the way, since the temperature of 300 K 
permits the thermal excitation of J = 4 of HD, then the 
rotational emission of HD (J=4-3) becomes possible. This 
emission occurs since the gas temperature can reach 300 K 
because of the H2 cooling balancing the adiabatic heating 
of the gravitationally contracting core. This efficiency of H2 
cooling is realized only when the three-body reaction among 
H atoms occurs. Thus, the second peak is possible owing to 
the three-body reaction to form H2. 



7 SUMMARY 

One of the main goals of cosmology is to find the first gener- 
ation of stars. When they form, strong H2 emission and weak 
HD emission is expected as a double peak feature. The weaJj 
HD emission is important since it distinguishes between the 
H2 emission from the primordial molecular clouds and that 
from the primordial supernova shell (Ciardi & Ferrara 2001). 
Wc have examined the observational feasibility of the detec- 
tion of the double peak emission. According to our analysis, 
the double peak feature is marginally detectable by ALMA. 
However, the transmissivity of air for the expected typical 
emission is low for the ALMA project as pointed out in 
our previous paper. If future telescopes are able to detect 
the double peak feature emission from primordial molecular 
cloud cores, then either primordial cores and the first stars 
have formed at different redshift from z = 19. 
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APPENDIX 

We re- formulate the energy loss due to lino emission of H2. 
The details are found in Kamaya & Silk (2001). Line emis- 
sion of H2 occurs due to the changes among rotation and 
vibration states. Wo can basically use the formulation of 
HoUenbach & McKee (1979) for rotational and vibrational 
emission of H2. Adopting their notation, we get: 



9.5 X 10" 



+3.0 X lO'^^exp 



-exp 
0.51 



'Oisy 



erg s 



(Al) 



-3 -1 
erg cm s , 



then, wo estimate the cooling rate as 

A(rOt) = nH2ir(l-|-CHr)"^+J'iH2ir(l-|-CH2i 

(A2) 

where (ht = nHcd(rot)/nH, Cngr = nH2cd(rot)/nH2, 
nHcd(rot) = Aj/7j, nH2cd(rot) = Aj/^^'^, and Aj is the 
Einstein A value for the J to J — 2 transition; 7j is the col- 
lisional de-excitation rate coefficient due to neutral hydro- 
gen; and is that due to molecular hydrogen. The first 
term of Lr denotes the cooling coefficient due to the higher 
rotation level (J > 2) and the second one due to J = 2 — » 
transition. The vibrational levels of both terms are set to be 
zero. 

For the vibrational transitions; 



Lv = 6.7xl0""exp 



-t-LGxlO-'^^exp 



(A3) 



then, we get the cooling rate as being 

A(vib) = nH2-I'v(l+CHv)~^+nH2-i'v(l+CH2v)~^ ergcm~^s~^ 

(A4) 

where Chv = nHcd(vib)/nH, and Ch2v = nH2cd(vib)/nH2- 
Here, nHcd(vib) = Aij/^^j, and nH2cd(vib) = Aij/^^^ where 
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Aij is the Einstein A value for the i to j transition; 7^ is 
the coUisional de-excitation rate coefficient due to neutral 
hydrogen; and 'yf^ is that due to molecular hydrogen. In 
our formula, only the levels of v = 0, 1 and 2 are considered. 
This is sufficient since the temperature is lower than 2000 
K. The first term of Lv is a cooling coefficient oi 5v = 1 
and the second term is that of 5v = 2. The second term 
has effect on only the central region of our structural model, 
then it has no significant contribution to our conclusion. 
Combining Eq.(A2) and Eq.(A4), wo obtain the total cooling 
rate as A"''" = A(rot)-|-A(vib) erg cm s ^ in the optically 
thin regime. When we need a cooling function which can be 
used in the optically thick regime, A"^'" is multiplied by the 
escape probability like the other cooling function in the main 
text. 

It may be better to comment on the continuum absorp- 
tion. The effect of the continuum absorption below 2000 K 
would be multiplied by exp(— rcont), in which 



'^cont 



^ p(r)Aj 



1 

p(r) 



1 

Po 



4.1 -f- - - Ts-'^ '" + 0.012p"-^' (r)T| 



(A5) 



according to the estimate of Lenzuni, Chernoff, & Salpeter 

(1991) who obtain a fitting formula for the Rossland mean 
opacity in a zero-metallicity gas. Here, Aj is the Jeans length 
and po is 0.8 g cm~^. Their fitting formula is reasonable if 
we consider the temperature range T > 1000 K. Our lowest 
temperature of the collapsing core is about 50 K. Then, their 
formula may not be appropriate, while it gives an sufficiently 
upper limit if we adopt Tcont of 1000 K instead of really 
having Tcont below 1000 K. For all of our estimates, Tcont is 
much smaller than unity. Then, we can neglect the effect of 
continuum absorption. 

Wc summarise the parameters in our calculation for a 
rotational transition with v = 0; A2,o = 2.94 x 10~^^ sec~^; 
and Aj,j-2 = 5^2,o/162 x J(J - 1)(2J - l)*/(2j + 1) sec"^ 
are considered. For a vibrational transition, Au) — 8.3 x 
10"^ sec"^; A21 = 1.1 X 10"^ scc"^; and A2Q = 4.1 x 10"'^ 
sec~^ are considered. In the main text, we find that the 
rotational line cooling of H2 is dominated by J = 2 — 
transition (« = 0). Then, we regard that the second term of 
Eq.(Al) is important. This means our estimate for ,7 = 2 — 
transition is very reliable, while the other estimate has some 
uncertainty. 

Finally, we comment on the uncertainty of the cooling 
function owing to rot-vibrational transitions of H2. The for- 
mula for the cooling function of H2 is examined by Martin 
et al. (1996), Forrey et al. (1997), GaUi & Palla (1998), and 
Puller &: Couchman (2000). Fuller & Couchman especially 
stress that there is uncertainty in the H2 cooling function 
because of the difficulty in calculating the interaction po- 
tential at low temperatures. Then, different choices for ro- 
tational and vibrational H-H2 rate coefficients will produce 
differences in the cooling function. Fortunately, we can con- 
sider our cooling function to be applicable since the cooling 
rate via H2 transitions balances the release of gravitational 
potential energy consistently as shown in our first paper 
(Kamaya & Silk 2001). 



REFERENCES 

Abel T., Bryan G. L., Norman M. L., 2000, ApJ, 540, 39 
Abgrall H., RoucflF R., Roueff E., 1992, A&AS, 50, 505 
Bougleux E., Galli D., 1997, MNRAS, 288, 638 
Bromm V., Coppi P. S., Larson R. B., 1999, ApJ, 527, L5 
Bromm V., Coppi P. S., Larson R. B., 2001, ApJ, 564, 23 
Carr B. J., Bond J. R., Arnett W. D., 1984, ApJ, 277, 445 
Ciardi B., Ferrara A., 2001, MNRAS, 324, 648 
Combes F., Pfenniger D., 1997, A&A, 327, 453 
Couchman H., Rees M., 1986, MNRAS, 221, 445 
Ferrara A., 1998, ApJ, 499, L17 

Forrey R. C, BalaJjrishnan N., Dalgarno A., Lepp S., 1997, ApJ, 

489, 100 

Fuller T. M., Couchman H. M. P., 2000, ApJ, 544, 6 

Galli D., Palla F., 1998, A&A, 335, 403 

Gncdin N. Y., Ostrlker J. P., 1997, ApJ, 486, 581 

Haiman Z., Thoul A. A., Loeb A., 1996, ApJ, 464, 523 

Hirashita H., Hunt L. K., Ferrara A., 2002, MNRAS, 330, L19 

Hollenbach D., McKee C. F., 1979, ApJS, 41, 555 

Kamaya H., Silk J., 2001, MNRAS, 332, 251 

Lenzuni P., Chcrnoflf D. F., Salpctcr E. E., 1991, ApJS, 76, 759 

Lcpp D., ShuU J. M., 1984, ApJ, 270, 578 

Low C, Lyndcn-Bell D., 1976, MNRAS, 176, 367 

Martin P. G., Schwarz D. H., Mandy M. E., 1996, ApJ, 461, 265 

Matsuda T., Sato H., Takcda H., 1969, Prog.Thcor.Phys., 41, 840 

Nishi R., Susa H., Uehara H., Yamada M., Omukai K., 1998, 

Prog.Theor.Phys., 100, 881 
Omukai K., Nishi R., 1998, ApJ, 508, 141 
Omukai K., 2000, ApJ 534, 809 

Palla F., Salpctcr E. E., Stabler W., 1983, ApJ, 271, 632 

Peebles P. J., Dicke R. H., 1968, ApJ, 154, 891 

Ripamonti, E., Haardt, F., Ferrara, A., Colpi, M., 2002, MNRAS, 

334, 401 

Saslaw W. C, Zipoy D., 1968, Nature, 216, 976 

Schaefer J., 1990, ApJS, 85, 1101 

Shchekinov lu., 1991, Ap&SS, 175, 57 

Shchekinov Yu., 1986, SvAL, 12, 211 

Silk J., 1977, ApJ, 214, 152 

Silk J., 1983, MNRAS, 205, 705 

Silk J., 2000, PASP, 112, 1003 

Susa H., Umcmura M., 2000, MNRAS, 316, L17 

Suto Y., Silk J., 1988, ApJ, 326, 527 

Takeuchi T. T., Kawabc R., Kohno K., Nakanishi K., Ishii T. T., 
Hirashita H., Yoshikawa K., 2001, PASP, 113, 586 

Toleman, R. C, 1966, Relativity, thermodynamics and cosmology, 
(Oxford; Clarendon Press), section 10 

Trefler M., Gush H.P., 1968, Phys.Rev.Let., 20, 703 

Uehara H., Inutsuka S., 2000, ApJ, 531, L91 



8 H.Kamaya and J. Silk 



Table 1. Expected Emission Lines 





H2 

( J = 2 - 0) 


HD 

(J = 4-3) 


LiH 

( J = 17 - 16) 






Di9 (lOio pc) 


u (IQl^ Hz; 2 = 0) 


1.06 


1.06 


0.74 








1/ (10l2 Hz; 2 = 19) 


0.53 


0.53 


0.37 








A (10"-^ mm; 2 = 0) 


28.29 


28.29 


40.30 








A (10~3 mm; z = 19) 


565.8 


565.8 


806.1 








(lO^^ erg sec-l Hz"!) 


40.1 


2.01 


0.001 








U (10-8 Jy; z = 19) 


43.18 


2.16 


0.0003 


1.0 


0.0 


0.62 


{10-« Jy; z = 19) 


16.60 


0.83 


0.0001 


0.3 


0.7 


1.00 


(10-« Jy; z = 19) 


7.58 


0.37 


0.00004 


0.1 


0.9 


1.48 
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Struclual Model 




log density (cm-^ 



Figure 1. Density-H2 fraction relation (solid line) and Density- 
Temperature relation (dashed line) are depicted. Using Eq.(8), we 
translate them to Radius-Density and -Temperature relations, 
respectively. 



Cooling Rate 




Figure 3. Heating rate; compressional (solid), reactive (dashed), 
and total (dotted). 



Cooling-Heating Ratio 



p 




log density (cm"^ 



Figure 4. Ratio of cooling to heating. As found clearly, our 
model is reasonable almost a few factors. This supports our model 
structure of a molecular cloud core. 



log density (cm"^) 

Figure 2. Cooling rate; HD (solid), H2 (short dashed), LiH (long 
dashed), and total (dotted). Below ~ 10* cm~^, HD is dominant, 
while H2 is important above the density. LiH is not so significant 
because of the lines of LiH is optically thick. 



Schematic View of Line Emission 
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Figure 5. Schematic view of the double peaJc emission of H2 
and HD. The diflference between the two lines is about 5.0x10^" 
Hz. This is resolved by ALMA project. 



